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ABSTRACT 

During the period of reionization the Universe was filled with a cosmological back- 
ground of ionizing radiation. By that time a significant fraction of the cosmic gas had 
already been incorporated into collapsed galactic halos with virial temperatures <^ 10 4 K 
that were unable to cool efficiently. We show that photoionization of this gas by the 
fresh cosmic UV background boiled the gas out of the gravitational potential wells of its 
host halos. We calculate the photoionization heating of gas inside spherically symmetric 
dark matter halos, and assume that gas which is heated above its virial temperature 
is expelled. In popular Cold Dark Matter models, the Press-Schechter halo abundance 
implies that ~ 50-90% of the collapsed gas was evaporated at reionization. The gas 
originated from halos below a threshold circular velocity of ~ 10-15 km s _1 . The re- 
sulting outflows from the dwarf galaxy population at redshifts z = 5-10 affected the 
metallicity, thermal and hydrodynamic state of the surrounding intergalactic medium. 
Our results suggest that stellar systems with a velocity dispersion <^ 10 km s _1 , such 
as globular clusters or the dwarf spheroidal galaxies of the Local Group, did not form 
directly through cosmological collapse at high redshifts. 



Subject headings: cosmology:theory — galaxiesdormation — galaxies: halos — radiative 
transfer 



1. Introduction 

The formation of galaxies is one of the most important, yet unsolved, problems in cosmology. 
The properties of galactic dark matter halos are shaped by gravity alone, and have been rigorously 
parameterized in hierarchical Cold Dark Matter (CDM) cosmologies (e.g., Navarro, Frenk, & White 
1997). However, the complex processes involving gas dynamics, chemistry and ionization, and 
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cooling and heating, which are responsible for the formation of stars from the baryons inside these 
halos, have still not been fully explored theoretically. 

Recent theoretical investigations of early structure formation in CDM models have led to a 
plausible picture of how the formation of the first cosmic structures leads to reionization of the 
intergalactic medium (IGM). The bottom-up hierarchy of CDM cosmologies implies that the first 
gaseous objects to form in the Universe have a low-mass, just above the cosmological Jeans mass of 
~ 1O 4 M (see, e.g., Haiman, Thoul, & Loeb 1996, and references therein). The virial temperature 
of these gas clouds is only a few hundred K, and so their metal-poor primordial gas can cool only 
due to the formation of molecular hydrogen, H2. However, H2 molecules are fragile, and were easily 
photo-dissociated throughout the Universe by trace amounts of starlight (Stecher & Williams 1967; 
Haiman, Rees, & Loeb 1996) that were well below the level required for complete reionization of 
the IGM. Following the prompt destruction of their molecular hydrogen, the early low-mass objects 
maintained virialized gaseous halos that were unable to cool or fragment into stars. Most of the 
stars responsible for the reionization of the Universe formed in more massive galaxies, with virial 
temperatures T v ; r ^ 10 4 K, where cooling due to atomic transitions was possible. The corresponding 
mass of these objects at z ~ 10 was ~ 10 s M Q , typical of dwarf galaxies. 

The lack of a Gunn-Peterson trough and the detection of Lya emission lines from sources out 
to redshifts z = 5.6 (Weymann et al. 1998; Dey et al. 1998; Spinrad et al. 1998; Hu, Cowie, &z 
McMahon 1998) demonstrates that reionization due to the first generation of sources must have 
occurred at yet higher redshifts; otherwise, the damping wing of Lya absorption by the neutral 
IGM would have eliminated the Lya line in the observed spectrum of these sources (Miralda- 
Escude 1998). Popular CDM models predict that most of the intergalactic hydrogen was ionized 
at a redshift 8 <; z <; 15 (Gnedin & Ostriker 1997; Haiman & Loeb 1998a,c). The end of the 
reionization phase transition resulted in the emergence of an intense UV background that filled the 
Universe and heated the IGM to temperatures of ~ 1-2 x 10 4 K (Haiman & Loeb 1998b; Miralda- 
Escude, Haehnelt, & Rees 1998). After ionizing the rarefied IGM in the voids and filaments on 
large scales, the cosmic UV background penetrated the denser regions associated with the virialized 
gaseous halos of the first generation of objects. Since a major fraction of the collapsed gas had been 
incorporated by that time into halos with a virial temperature <^ 10 4 K, photoionization heating by 
the cosmic UV background could have evaporated much of this gas back into the IGM. No such 
feedback was possible at earlier times, since the formation of internal UV sources was suppressed 
by the lack of efficient cooling inside most of these objects. 

The gas reservoir of dwarf galaxies with virial temperatures ^ 10 4 K (or equivalently a ID 
velocity dispersion ^ 10 km s" 1 ) could not be immediately replenished. The suppression of dwarf 
galaxy formation at z > 2 has been investigated both analytically (Rees 1986; Efstathiou 1992) and 
with numerical simulations (Thoul & Weinberg 1996; Quinn, Katz, k, Efstathiou 1996; Weinberg, 
Hernquist, & Katz 1997; Navarro & Steinmetz 1997). The dwarf galaxies which were prevented 
from forming after reionization could have eventually collected gas at z = 1-2, when the UV 
background flux declined sufficiently (Babul & Rees 1992; Kepner, Babul, & Spergel 1997). The 
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reverse process during the much earlier reionization epoch has not been addressed in the literature. 
(However, note that the photo-evaporation of gaseous halos was considered by Bond, Szalay, & Silk 
(1988) as a model for Lya absorbers at lower redshifts z ~ 4.) 

In this paper we focus on the reverse process by which gas that had already settled into 
virialized halos by the time of reionization was evaporated back into the IGM due to the cosmic 
UV background which emerged first at that epoch. The basic ingredients of our model are presented 
in §2. In order to ascertain the importance of a self-shielded gas core, we include a realistic, centrally 
concentrated dark halo profile and also incorporate radiative transfer. Generally we find that self- 
shielding has a small effect on the total amount of evaporated gas, since only a minor fraction of 
the gas halo is contained within the central core. Our numerical results are described in §3. In 
particular, we show the conditions in the highest mass halo which can be disrupted at reionization. 
We also use the Press-Schechter (1974) prescription for halo abundance to calculate the fraction 
of gas in the Universe which undergoes the process of photo-evaporation. Our versatile semi- 
analytic approach has the advantage of being able to yield the dependence of the results on a wide 
range of reionization histories and cosmological parameters. Clearly, the final state of the gas halo 
depends on its dynamical evolution during its photo-evaporation. We adopt a rough criterion for 
the evaporation of gas based on its initial interaction with the ionizing background. The precision 
of our results could be tested in specific cases by future numerical simulations. In §4 we discuss the 
potential implications of our results for the state of the IGM and for the early history of low-mass 
galaxies in the local Universe. Finally, we summarize our main conclusions in §5. 



A Model for Halos at Reionization 



We consider gas situated in a virialized dark matter halo. We adopt the prescription for ob- 
taining the density profiles of dark matter halos at various redshifts from the Appendix of Navarro, 
Frenk, & White (1997, hereafter NFW), modified to include the variation of the collapse overdensity 
A c . Thus, a halo of mass M at redshift z is characterized by a virial radius, 
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The density profile of the halo is given by 

( r ) = ^1 (l + z) 3 - 
P{ ' 8irG [ ^ ' n(z) cx(l + cxf ' 

where x = r/r v - a and c depends on <5 C for a given mass M. We include the dependence of halo 
profiles on fi and Q.\, the current contributions to 17 from non-relativistic matter and a cosmological 
constant, respectively (see Appendix A for complete details). 
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Although the NFW profile provides a good approximation to halo profiles, there are indications 
that halos may actually develop a core (e.g., Burkert 1995; Kravtsov et al. 1998; see, however, Moore 
et al. 1998). In order to examine the sensitivity of the results to model assumptions, we consider 
several different gas and dark matter profiles, keeping the total gas fraction in the halo equal to 
the cosmological baryon fraction. The simplest case we consider is an equal NFW profile for the 
gas and the dark matter. In order to include a core, instead of the NFW profile of equation (|3j) we 
also consider the density profile of the form fit by Burkert (1995) to dwarf galaxies, 

P(r) = 8^G (1 + Z)3 n(z) (1 + bx) [1 + [bxf] (4) 

where b is the inverse core radius, and we set S c by requiring the mean overdensity to equal the 
appropriate value, A c , in each cosmology (see Appendix A). We also consider two cases where the 
dark matter follows an NFW profile but the gas is in hydrostatic equilibrium with its density profile 
determined by its temperature distribution. In one case, we assume the gas is isothermal at the 
halo virial temperature, given by 

rp _VV C 2 _ ocinn M ( r vh \ 2 n A c {l + z\ 3 
vir " 2k B ~ 6bLW 0.6m p U^kpJ Q(z)20o{ 10 ) ' W 

where \x is the mean molecular weight as determined by ionization equilibrium, and m p is the 
proton mass. The spherical collapse simulations of Haiman, Thoul, & Loeb (1996) find a post- 
shock gas temperature of roughly twice the value given by equation (^), so we also compare with 
the result of setting T = 2 T v i r . In the second case, we let the gas cool for a time equal to the 
Hubble time at the redshift of interest, z. Gas above 10 4 K cools rapidly due to atomic cooling 
until it reaches a temperature near 10 4 K, where the cooling time rapidly diverges. In this case, 
hydrostatic equilibrium yields a highly compact gas cloud when the halo virial temperature is 
greater than 10 4 K. In reality, of course, a fraction of the gas may fragment and form stars in these 
halos. However, this caveat hardly affects our results since only a small fraction of the gas which 
evaporates is contained in halos with T vir > 10 4 K. Throughout most of our subsequent discussion 
we consider the simple case of identical NFW profiles for both the dark matter and the gas, unless 
indicated otherwise. 

We assume a helium mass fraction of Y = 0.24, and include it in the calculation of the ionization 
equilibrium state of the gas as well as its cooling and heating (see, e.g., Katz, Weinberg, & Hernquist 
1996). We adopt the various reaction and cooling rates from the literature, including the rates for 
collisional excitation and dielectronic recombination from Black (1981); the recombination rates 
from Verner & Ferland (1996), and the recombination cooling rates from Ferland et al. (1992) with 
a fitting formula by Miralda-Escude (1998, private communication). Collisional ionization rates are 
adopted from Voronov (1997), with the corresponding cooling rate for each atomic species given by 
its ionization rate multiplied by its ionization potential. We also include cooling by Bremsstrahlung 
emission with a Gaunt factor from Spitzer & Hart (1971), and by Compton scattering off the 
microwave background (e.g., Shapiro &: Kang 1987). 
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In assessing the effect of reionization, we assume for simplicity a sudden turn-on of an external 
radiation field with a specific intensity per unit frequency, v, 



where vl is the Lyman limit frequency. Our treatment of the response of the cloud to this radiation, 
as outlined below, is not expected to yield different results with a more gradual increase of the 
intensity with cosmic time. The external intensity is responsible for the reionization of the 

IGM, and so we normalize it to have a fixed number of ionizing photons per baryon in the Universe. 
We define the ionizing photon density as 



where the photoionization efficiency is weighted by the photoionization cross section of HI, ghi{v), 
above the Lyman limit. The mean baryon number density is 



Throughout the paper we refer to proper densities rather than comoving densities. As our standard 
case we assume a post-reionization ratio of n 7 /n^ = 1, but we also consider the effect of setting 
n 7 /nb = 0.1. For example, a = 1.8 and n 7 /nj, = 1 yield I^i = 1.0 at z = 3 and I%i = 3.5 at 
z = 5, close to the values required to satisfy the Gunn-Peterson constraint at these redshifts (see, 
e.g., Efstathiou 1992). Note that n^/rib ^ 1 is required for the initial ionization of the gas in the 
Universe (although this ratio may decline after reionization). 

We assume that the above uniform UV background illuminates the outer surface of the gas 
cloud, located at the virial radius r v i r , and penetrates from there into the cloud. The radiation 
photoionizes and heats the gas at each radius to its equilibrium temperature, determined by equat- 
ing the heating and cooling rates. The latter assumption is justified by the fact that both the 
recombination time and the heating time are initially shorter than the dynamical time throughout 
the halo. At the outskirts of the halo the dynamics may start to change before the gas can be 
heated up to its equilibrium temperature, but this simply means that the gas starts expanding out 
of the halo during the process of photoheating. This outflow should not alter the overall fraction 
of evaporated gas. 

The process of reionization is expected to be highly non-uniform due to the clustering of the 
ionizing sources and the dumpiness of the IGM. As time progresses, the HII regions around the ion- 
izing sources overlap, and each halo is exposed to ionizing radiation from an ever increasing number 
of sources. While the external ionizing radiation may at first be dominated by a small number of 
sources, it quickly becomes more isotropic as its intensity builds up with time (e.g., Haiman & Loeb 
1998a,b; Miralda-Escude, Haehnelt, & Rees 1998). The evolution of this process depends on the 
characteristic clustering scale of ionizing sources and their correlation with the inhomogeneities of 
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the IGM. In particular, the process takes more time if the sources are typically embedded in dense 
regions of the neutral IGM which need to be ionized first before their radiation shines on the rest 
of the IGM. However, in our analysis we do not need to consider these complications since the total 
fraction of evaporated gas in bound halos depends primarily on the maximum intensity achieved 
at the end of the reionization epoch. 

In computing the effect of the background radiation, we include self-shielding of the gas which 
is important at the high densities obtained in the core of high redshift halos. For this purpose, 
we include radiative transfer through the halo gas and photoionization by the resulting anisotropic 
radiation field in the calculation of the ionization equilibrium. We also include the fact that the 
ionizing spectrum becomes harder at inner radii, since the outer gas layers preferentially block 
photons with energies just above the Lyman limit. We neglect self-shielding due to helium atoms. 
Appendix B summarizes our simplified treatment of the radiative transfer equations. 

Once the gas is heated throughout the halo, some fraction of it acquires a sufficiently high 
temperature that it becomes unbound. This gas expands due to the resulting pressure gradient and 
eventually evaporates back to the IGM. The pressure gradient force (per unit volume) fceV(Tp//x) 
competes with the gravitational force of pGM/r 2 . Due to the density gradient, the ratio between 
the pressure force and the gravitational force is roughly the ratio between the thermal energy ~ ksT 
and the gravitational binding energy ~ fiGM/r (which is ~ ksT^ at r v - ir ) per particle. Thus, if the 
kinetic energy exceeds the potential energy (or roughly if T > T v i r ), the repulsive pressure gradient 
force exceeds the attractive gravitational force and expels the gas on a dynamical time (or faster 
for halos with T > T vir ). 

We compare the thermal and gravitational energy (both of which are functions of radius) as 
a benchmark for deciding which gas shells are expelled from each halo. Note that infall of fresh 
IGM gas into the halo is also suppressed due to its excessive gas pressure, produced by the same 
photo-ionization heating process. 

This situation stands in contrast to feedback due to supernovae, which depends on the effi- 
ciency of converting the mechanical energy of the supernovae into thermal energy of the halo gas. 
The ability of supernovae to disrupt their host dwarf galaxies has been explored in a number of 
theoretical papers (e.g., Larson 1974; Dekel & Silk 1986; Vader 1986, 1987). However, numerical 
simulations (Mac-Low & Ferrara 1998) find that supernovae produce a hole in the gas distribution 
through which they expel the shock-heated gas, leaving most of the cooler gas still bound. In 
the case of reionization, on the other hand, energy is imparted to the gas directly by the ionizing 
photons. A halo for which a large fraction of the gas is unbound by reionization is thus prevented 
from further collapse and star formation. 

When the gas in each halo is initially ionized, an ionization shock front may be generated (cf. 
the discussion of Lya absorbers by Donahue & Shull 1987). The dynamics of such a shock front 
have been investigated in the context of the interstellar medium by Bertoldi & McKee (1990) and 
Bertoldi (1989). Their results imply that the dynamics of gas in a halo are not significantly affected 
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by the shock front unless the thermal energy of the ionized gas is greater than its gravitational 
potential energy. Furthermore, since gas in a halo is heated to the virial temperature even before 
reionization, the shock is weaker when the gas is ionized than a typical shock in the interstellar 
medium. Also, as noted above, the ionizing radiation reaching a given halo builds up in intensity 
over a considerable period of time. Thus, we do not expect the ionization shock associated with the 
first encounter of ionizing radiation to have a large effect on the eventual fate of gas in the halo. 

3. Results 

We assume the most popular cosmology to date (Garnavich et al. 1998) with f^o = 0.3 and 
f&A = 0.7. We illustrate the effects of cosmological parameters by displaying the results also for 
Qq = 1, and for = 0.3 and J7a = 0. The models all assume f^/i 2 = 0.02 and a Hubble constant 
h = 0.5 if Qq = 1 and h = 0.7 otherwise (where Hq = 100 h km s~ 1 Mpc _1 ). 

Figure 1 shows the temperature of the gas versus its baryonic overdensity A& relative to the 
cosmic average (cf. Efstathiou 1992). The curves are for z = 8 and assume £Iq = 0.3 and £l\ = 0.7. 
We include intergalactic radiation with a flux given by equation (||) for a = 1.8 and n^/ri}, = 1. 
The dotted curve shows tjj = t coo i with no radiation field, where tu is the age of the Universe, 
approximately equal to 6.5 x 10 9 /i _1 (l + z)~ 3 / 2 Sl years at high redshift. This curve indicates 
the temperature to which gas has time to cool through atomic transitions before reionization. This 
temperature is always near T = 10 4 K since below this temperature the gas becomes mostly neutral 
and the cooling time is very long. It is likely that only atomic cooling is relevant before reionization 
since molecular hydrogen is easily destroyed by even a weak ionizing background (Haiman, Rees, 
& Loeb 1996). The solid curve shows the equilibrium temperature for which the heating time 
theat due to a UV radiation field equals the cooling time t coo i. The decrease in the temperature at 
Aft < 10 is due to the increased importance of Compton cooling, which is proportional to the gas 
density rather than its square. At a given density, gas is heated at reionization to the temperature 
indicated by the solid curve, unless the net cooling or heating time is too long. The dashed curves 
show the temperature where the net cooling or heating time equals tu. By definition, points on 
the solid curve have an infinite net cooling or heating time, but there is also a substantial regime 
at low A;, where the net cooling or heating time is greater than tu. However, this regime has only 
a minor effect on halos, since the mean overdensity inside the virial radius of a halo is of order 200. 
On the other hand, if gas leaves the halo and expands it quickly enters the regime where it cannot 
reach thermal equilibrium. 

Figure 2 presents an example for the structure of a halo with an initial total mass of M = 
3 x 10 7 Mq at z = 8. We assume the same cosmological parameters as in Figure 1. The bottom 
plot shows the baryon overdensity A& versus r/r v ; r , and reflects our assumption of identical NFW 
profiles for both the dark matter and the baryons. The middle plot shows the neutral hydrogen 
fraction versus r/r v i r , and the top plot shows the ratio of thermal energy per particle (TE = §fc fl T) 
to potential energy per particle (PE = fj,\<f)(r)\, where 0(r) is the gravitational potential) versus 
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r/r v ; r . The dashed curves assume an optically thin halo, while the solid curves include radiative 
transfer and self-shielding. The self-shielded neutral core is apparent from the solid curves, but 
since the point where TE/PE = 1 occurs outside this core, the overall unbound fraction does not 
depend strongly on the radiative transfer in this case. Its value is 67% assuming an optically-thin 
halo, and 64% when radiative transfer is included and only a fraction of the external photons make 
their way inside. Even when the opacity at the Lyman limit is large, some ionizing radiation still 
reaches the central parts of the halo because, (i) the opacity drops quickly above the Lyman limit, 
and (ii) the heated gas radiates ionizing photons inwards. 

Figure 3 shows the unbound gas fraction after reionization as a function of the total halo mass. 
We assume £Iq = 0.3, £1\ = 0.7, and n^/rib = 1. The three pairs of curves shown consist of a solid 
line (which includes radiative transfer) and a dashed line (which assumes an optically thin halo). 
From right to left, the first pair is for a = 1.8 and z = 8, the second is for a = 5 and z = 8, 
and the third is for a = 1.8 and z = 20. In each case the self-shielded core lowers the unbound 
fraction when we include radiative transfer (solid vs dashed lines), particularly when the unbound 
fraction is sufficiently large that it includes part of the core itself. High energy photons above the 
Lyman limit penetrate deep into the halo and heat the gas efficiently. Therefore, a steepening of the 
spectral slope from a = 1.8 to a = 5 decreases the temperature throughout the halo and lowers the 
unbound gas fraction. This is only partially compensated for by our UV flux normalization, which 
increases I21 with increasing a so as to get the same density of ionizing photons in equation (0). 
Increasing the reionization redshift from z = 8 to z = 20 increases the binding energy of the gas, 
because the high redshift halos are denser. Although the corresponding increase of I21 with redshift 
(at a fixed n 7 /n&) counteracts this change, the fraction of expelled gas is still reduced due to the 
deeper potential wells of higher redshift halos. 

From plots similar to those shown in Figure 3, we find the total halo mass at which the unbound 
gas fraction is 50%. We henceforth refer to this mass as the 50% mass. Figure 4 plots this mass 
as a function of the reionization redshift for different spectra and cosmological models. The solid 
line assumes a = 1.8 and the dotted line a = 5, both for Qq = 0.3 and = 0.7. The other 
lines assume a = 1.8 but different cosmologies. The short-dashed line assumes = 0.3, Qa = 
and the long-dashed line assumes = 1. All assume n^/rib = 1. Gas becomes unbound when its 
thermal energy equals its potential binding energy. The thermal energy depends on temperature, 
but the equilibrium temperature does not change much with redshift since we increase the UV flux 
normalization by the same (1 + z) 3 factor as the mean baryonic density. With this prescription for 
the UV flux, the 50% mass occurs at a value of the circular velocity which is roughly constant with 
redshift. Thus for each curve, the change in mass with redshift is mostly due to the change in the 
characteristic halo density, which affects the relation between circular velocity and mass. 

The cosmological parameters have only a modest effect on the 50% mass, and change it by up 
to 35% at a given redshift. Lowering reduces the characteristic density of a halo of given mass, 
and so a higher mass is required in order to keep the gas bound. Adding a cosmological constant 



reduces the density further through A c [see equations ( |10[) and (11)]. For the three curves with 
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a = 1.8, the circular velocity of the 50% mass equals 13 km s 1 at all redshifts, up to variations of 
a few percent. 

The spectral shape of the ionizing flux affects modestly the threshold circular velocity cor- 
responding to the 50% mass, because assuming a steeper spectrum (i.e. with a larger a) reduces 
the gas temperature and thus requires a shallower potential to keep the gas bound. A higher flux 
normalization has the opposite effect of increasing the threshold circular velocity. The left panel 
of Figure 5 shows the variation of circular velocity with spectral shape, for two normalizations 
(n 7 /n& = 1 and n^/n/j = 0.1 for the solid and dashed curves, respectively). The right panel shows 
the complementary case of varying the spectral normalization, using two values for the spectral 
slope (a = 1.8 and a = 5 for the solid and dashed curves, respectively). All curves assume an 
Qo = 0.3, = 0.7 cosmology. 

Obviously, 50% is a fairly arbitrary choice for the unbound gas fraction at which halos evapo- 
rate. Figure 3 shows that for a given halo, the unbound gas fraction changes from 10% to 90% over 
a factor of ~ 60 in mass, or a factor of ~ 4 in velocity dispersion. When 50% of the gas is unbound, 
however, the rest of the gas is also substantially heated, and we expect the process of collapse and 
fragmentation to be inhibited. In the extreme case where the gas expands until a steady state is 
achieved where it is pressure confined by the IGM, less than 10% of the original gas is left inside the 
virial radius. However, continued infall of dark matter should limit the expansion. Numerical sim- 
ulations may be used to define more precisely the point at which gas halos are disrupted. Clearly, 
photo-evaporation affects even halos with masses well above the 50% mass, although these halos 
do not completely evaporate. Note that it is also clear from Figure 3 that not including radiative 
transfer would have only a minor effect on the value of the 50% mass (typically ~ 5%). 

Given the values of the unbound gas fraction in halos of different masses, we can integrate to 
find the total gas fraction in the Universe which becomes unbound at reionization. This calculation 
requires the abundance distribution of halos, which is readily provided by the Press-Schechter mass 
function for CDM cosmologies (relevant expressions are given, e.g., in NFW). The high-mass cutoff 
in the integration is given by the lowest mass halo for which the unbound gas fraction is zero, since 
halos above this mass are not significantly affected by the UV radiation. The low-mass cutoff is 
given by the lowest mass halo in which gas has assembled by the reionization redshift. We adopt 
for this low-mass cutoff the linear Jeans mass, which we calculate following Peebles (1993, §6). The 
gas temperature in the Universe follows the cosmic microwave background temperature down to a 
redshift 1 + zt ~ 740($7b/i 2 ) 2//5 , at which the baryonic Jeans mass is 1.9 x 10 5 (ft fe /i 2 )~ 1/2 M o . After 
this redshift, the gas temperature goes down as (1 + z) 2 , so the baryon Jeans mass acquires a factor 
of [{l + z)/{l + z t )f/ 2 . Until now we have considered baryons only, but if we add dark matter then 
the mean density (or the corresponding gravitational force) is increased by fio/^bj which decreases 
the baryonic Jeans mass by — 3 / 2 . The corresponding total halo mass is f^o/^b times the 
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baryonic mass. Thus the Jeans cutoff before reionization corresponds to a total halo mass of 



We thus calculate the total fraction of gas in the Universe which is bound in pre-existing halos, 
and the fraction of this gas which then becomes unbound at reionization. In Figure 6 we show the 
fraction of the collapsed gas which evaporates as a function of the reionization redshift. The solid 
line assumes a. = 1.8, and the dotted line assumes a = 5, both for £Iq = 0.3, f^A = 0.7. The other 
lines assume a = 1.8, the short-dashed line with Qq = 0.3, £7a = and the long-dashed line with 
0,q = 1. All assume n^/n^ = 1 and a primordial n = 1 (scale invariant) power spectrum. In each 
case we normalized the CDM power spectrum to the present cluster abundance, a% = 0.5 Oq ' 5 
(see, e.g., Pen 1998), where as is the root-mean-square amplitude of mass fluctuations in spheres 
of radius 8 h^ 1 Mpc. The fraction of collapsed gas which is unbound is ~ 0.4-0.7 at z = 6 and 
it increases with redshift. This fraction clearly depends strongly on the halo abundance but is 
relatively insensitive to the spectral slope a of the ionizing radiation. In hierarchical models, the 
characteristic mass (and binding energy) of virialized halos is smaller at higher redshifts, and a 
larger fraction of the collapsed gas therefore escapes once it is photoheated. Among the three 
cosmological models, the characteristic mass at a given redshift is smallest for = 1 and largest 
for Q = 0.3, Q A = 0. 

In Figure 7 we show the total fraction of gas in the Universe which evaporates at reionization. 
The solid line assumes a = 1.8, and the dotted line assumes a = 5, both for Qq = 0.3, S7a = 0.7. 
The other lines assume a = 1.8, the short-dashed line with f^o = 0.3, f^A = and the long-dashed 
line with f^o = 1- All assume n^/rib = 1. For the different cosmologies, the total unbound fraction 
goes up to 20-25% if reionization occurs as late as z = 6-7; in this case a substantial fraction of 
the total gas in the Universe undergoes the process of expulsion from halos. However, this fraction 
typically decreases at higher redshifts. Although a higher fraction of the collapsed gas evaporates 
at higher z (see Figure 6), a smaller fraction of the gas in the Universe lies in halos in the first place. 
The latter effect dominates except for the open model up to z ~ 7. As is well known, the £1$ = 1 
model produces late structure formation, and indeed the collapsed fraction decreases rapidly with 
redshift in this cosmological model. The low Qq models approach the Oo = 1 behavior at high z, 
but this occurs faster for the flat model with a cosmological constant than for the open model with 
the same value of $7o- 

Changing the dark matter and gas profiles as discussed in §2 has a modest effect on the results. 
For example, with = 0.3, S7a = 0.7, and z = 8, and for our standard model where the gas and 
dark matter follow identical NFW profiles, the total unbound gas fraction is 19.8% and the halo 
mass which loses 50% of its baryons is 5.25 x 10 7 Mq. If we let the mass and the baryons follow 
the profile of equation (|j) the corresponding results are 20.0% and 5.31 x 10 7 Mq for b = 10 in 




0) 



This value agrees with the numerical spherical collapse calculations of Haiman, Thoul, & Loeb 
(1996). 
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equation (g) and 20.9% and 6.84 x 10 7 M Q for b = 5 ( i.e. a larger core). With an NFW mass profile 
but gas in hydrostatic equilibrium at the virial temperature, the unbound fraction is 19.2%, and 
the 50% mass is 4.33 x 1O 7 M . If we let the gas temperature be T = 2 T Y - ir , the unbound fraction 
is 22.0% and the 50% mass is 1.18 x 10 8 Mq. For clouds of gas which condense by cooling for a 
Hubble time, the unbound fraction is 18.2%, and the 50% mass is 3.38 x 1O 7 M . We conclude 
that centrally concentrated gas clouds are in general more effective at retaining their gas, but the 
effect on the overall unbound gas fraction in the Universe is modest, even for large variations in 
the profile. If we return to the NFW profile but adopt / = 0.01 instead of / = 0.5 in the NFW 
prescription for finding the collapse redshift (see Appendix A), we find an unbound fraction of 
20.3%, and a 50% mass of 6.06 x 10 7 Mq. Finally, lowering Qf, by a factor of 2 changes the unbound 
fraction to 19.0% and the 50% mass to 5.44 x 10 7 Mq. Our predictions appear to be robust against 
variations in the model parameters. 



4. Implications for the Intergalactic Medium and for Low Redshift Objects 

Our calculations show that a substantial fraction of gas in the Universe may lie in virialized 
halos before reionization, and that most of it evaporates out of the halos when it is photoionized and 
heated at reionization. The resulting outflows of gas from halos may have interesting implications 
for the subsequent evolution of structure in the IGM. We discuss some of these implications in this 
section. 

In the pre-reionization epoch, a fraction of the gas in the dense cores of halos may fragment 
and form stars. Some star formation is, of course, needed in order to produce the ionizing flux 
which leads to reionization. These population III stars produce the first metals in the Universe, 
and they may make a substantial contribution to the enrichment of the IGM. Numerical models 
by Mac-Low & Ferrara (1998) suggest that feedback from supernovae is very efficient at expelling 
metals from dwarf galaxies of total mass 3.5 x 10 s Mq, although it ejects only a small fraction of 
the interstellar medium in these hosts. Obviously, the metal expulsion efficiency depends on the 
presence of clumps in the supernova ejecta (Franco et al. 1993) and on the supernova rate - the 
latter depending on the unknown star formation rate and the initial mass function of stars at high 
redshifts. Reionization provides an alternative method for expelling metals efficiently out of dwarf 
galaxies by directly photoheating the gas in their halos, leading to its evaporation along with its 
metal content |] 

Gas which falls into halos and is expelled at reionization attains a different entropy than if 
it had stayed at the mean density of the Universe. Gas which collapses into a halo is at a high 



3 Note that we have assumed zero metallicity in calculating cooling. Even if some metals had already been mixed 
into the IGM, the metallicity of newly formed objects was likely too low to affect cooling since even at z ~ 3 the 
typical metallicity of the Lyman alpha forest has been observed to be < 0.01 solar (Songaila & Cowie 1996; Tytler 
et al. 1995). 
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overdensity when it is photoheated, and is therefore at a lower entropy than if it were heated to 
the same temperature at the mean cosmic density. However, the overall change in the entropy 
density of the IGM is small for two reasons. First, even at z = 6 only about 25% of the gas 
in the Universe undergoes evaporation. Second, the gas remains in ionization equilibrium and is 
photoheated during its initial expansion. For example, if z = 6, Qo = 0.3, S7a = 0.7, n 7 /n{, = 1, 
and a = 1.8, then the recombination time becomes longer than the dynamical time only when the 
gas expands down to an overdensity of 26, at which point its temperature is 22,400 K compared 
to an initial (non-equilibrium) temperature of 19,900 K for gas at the mean density. The resulting 
overall reduction in the entropy is the same as would be produced by reducing the temperature of 
the entire IGM by a factor of 1.6. This factor reduces to 1.4 if we increase z to 8 or increase a to 5. 
Note that Haehnelt & Steinmetz (1998) showed that differences in temperature by a factor of 3-4 
result in possibly observable differences in the Doppler parameter distribution of Lya absorption 
lines at redshifts 3-5. 

When the halos evaporate, recombinations in the gas could produce Lya lines or radiation 
from two-photon transitions to the ground state of hydrogen. However, a simple estimate shows 
that the resulting luminosity is too small for direct detection unless these halos are illuminated 
by an internal ionizing source. In an externally illuminated z = 6, 10 8 Mq halo our calculations 
imply a total of ~ 1 x 10 50 recombinations per second. Note that the number of recombinations 
is dominated by the high density core, and if we did not include self-shielding we would obtain an 
overestimate by a factor of ~ 15. If each recombination releases one or two photons with a total 
energy of 10.2 eV, then for S7o = 0.3 and S7a = 0.7 the observed flux is ~ 5 x 10 -20 erg s" 1 cm~ 2 . 
This flux is well below the sensitivity of the planned Next Generation Space Telescope, even if part 
of this flux is concentrated in a narrow line. 

The photoionization heating of the gaseous halos of dwarf galaxies resulted in outflows with 
a characteristic velocity of ~ 20-30 km s _1 . These outflows must have induced peculiar velocities 
of a comparable magnitude in the IGM surrounding these galaxies. The effect of the outflows 
on the velocity field and entropy of the IGM at z = 5-10 could in principle be searched for in 
the absorption spectra of high redshift sources, such as quasars. These small-scale fluctuations in 
velocity and the resulting temperature fluctuations have been seen in recent simulations by Bryan 
et al. (1998). However, the small halos responsible for these outflows were only barely resolved 
even in these high resolution simulations of a small volume. 

The evaporating galaxies could contribute to the high column density end of the Lya forest (cf. 
Bond, Szalay, Sz Silk 1988). For example, shortly after being photoionized, a z = 8, 5x 10 7 Mq halo 
has a neutral hydrogen column density of 2 x 10 16 cm -2 at an impact parameter of 0.5r v i r = 0.66 
kpc, 6xl0 17 cm" 2 at 0.25 r vir , and9xl0 20 cm" 2 (or 9xl0 18 cm if we do not include self-shielding) 
at 0.1r v i r (assuming Qq = 0.3, = 0.7, a = 1.8, and n^jn^ = 1). These column densities will 
decline as the gas expands out of the host galaxy. Abel & Mo (1998) have suggested that a 
large fraction of the Lyman limit systems at z ~ 3 may correspond to mini-halos that survived 
reionization. Remnant absorbers due to galactic outflows can be distinguished from large-scale 
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absorbers in the IGM by their compactness. Close lines of sight due to quasar pairs or gravitational 
lensed quasars (see, e.g., Crotts & Fang 1998; Petry, Impey, & Foltz 1998, and references therein) 
should probe different HI column densities in galactic outflow absorbers but similar column densities 
in the larger, more common absorbers. Follow-up observations with high spectroscopic resolution 
could reveal the velocity fields of these outflows. 

Although much of the gas in the Universe evaporated at reionization, the underlying dark 
matter halos continued to evolve through infall and merging, and the heated gas may have accu- 
mulated in these halos at lower redshifts. This latter process has been discussed by a number of 
authors, with an emphasis on the effect of reionization and the resulting heating of gas. Thoul & 
Weinberg (1996) found a reduction of ~ 50% in the collapsed gas mass due to heating, for a halo 
of V c = 50 km s _1 at z = 2, and a complete suppression of infall below V c = 30 km s _1 . The 
effect is thus substantial on halos with virial temperatures well above the gas temperature. Their 
interpretation is that pressure support delays turnaround substantially and slows the subsequent 
collapse. Indeed, as noted in §2, the ratio of the pressure force to the gravitational force on the gas 
is roughly equal to the ratio of its thermal energy to its potential energy. For a given enclosed mass, 
the potential energy of a shell of gas increases as its radius decreases. Before collapse, each gas shell 
expands with the Hubble flow until its expansion is halted and then reversed. Near turnaround, 
the gas is weakly bound and the pressure gradient may prevent collapse even for gas below the 
halo virial temperature. On the other hand, gas which is already well within the virial radius is 
tightly bound, which explains our lower value of V c ~ 13 km s~ x for halos which lose half their gas 
at reionization. 

Three dimensional numerical simulations (Quinn, Katz, & Efstathiou 1996; Weinberg, Hern- 
quist, k, Katz 1997; Navarro & Steinmetz 1997) have also explored the question of whether dwarf 
galaxies could re-form at z ^> 2. The heating by the UV background was found to suppress infall 
of gas into even larger halos (V c ~ 75 km s -1 ), depending on the redshift and on the ionizing 
radiation intensity. Navarro & Steinmetz (1997) noted that photoionization reduces the cooling 
efficiency of gas at low densities, which suppresses further the late infall at redshifts below 2. We 
note that these various simulations assume an isotropic ionizing radiation field, and do not calcu- 
late radiative transfer. Photoevaporation of a gas cloud has been calculated in a two dimensional 
simulation (Shapiro, Raga, & Mellema 1998), and methods are being developed for incorporating 
radiative transfer into three dimensional cosmological simulations (e.g., Abel, Norman, & Madau 
1999; Razoumov & Scott 1999). 

Our results have interesting implications for the fate of gas in low-mass halos. Gas evaporates 
at reionization from halos below V c ~ 13 km s _1 , or a velocity dispersion a ~ 10 km s _1 . A similar 
value of the velocity dispersion is also required to reach a virial temperature of 10 4 K, allowing 
atomic cooling and perhaps star formation before reionization. Thus, halos with a ^ 10 km s _1 
could have formed stars before reionization. They would have kept their gas after reionization, and 
could have had ongoing star formation subsequently. These halos were the likely sites of population 
III stars, and could have been the progenitors of dwarf galaxies in the local Universe (cf. Miralda- 
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Escude & Rees 1998). On the other hand, halos with a <^ 10 km s could not have cooled before 
reionization. Their warm gas was completely evaporated from them at reionization, and could not 
have returned to them until very low redshifts, possibly z ^ 1, so that their stellar population 
should be relatively young. 

It is interesting to compare these predictions to the properties of dwarf spheroidal galaxies in 
the Local Group which have low central velocity dispersions. At first sight this appears to be a 
difficult task. The dwarf galaxies vary greatly in their properties, with many showing evidence for 
multiple episodes of star formation as well as some very old stars (see the recent review by Mateo 
1998). Another obstacle is the low temporal resolution of age indicators for old stellar populations. 
For example, if Qo = 0.3 and = 0.7 then the age of the Universe is 43% of its present age at 
z = 1 and 31% at z = 1.5. Thus, stars that formed at these redshifts may already be ~ 10 Gyr old 
at present, and are difficult to distinguish from stars that formed at z > 5. 

Nevertheless, one of our robust predictions is that most early halos with a <^ 10 km s _1 could 
not have formed stars in the standard hierarchical scenario. Globular clusters belong to one class 
of objects with such a low velocity dispersion. Peebles & Dicke (1968) originally suggested that 
globular clusters may have formed at high redshifts, before their parent galaxies. However, in 
current cosmological models, most mass fluctuations on globular cluster scales were unable to cool 
effectively and fragment until z ~ 10, and were evaporated subsequently by reionization. We note 
that Fall & Rees (1985) proposed an alternative formation scenario for globular clusters involving 
a thermal instability inside galaxies with properties similar to those of the Milky Way. Globular 
clusters have also been observed to form in galaxy mergers (e.g., Miller et al. 1997). It is still 
possible that some of the very oldest and most metal poor globular clusters originated from z ^ 10, 
before the UV background had become strong enough to destroy the molecular hydrogen in them. 
However, primeval globular clusters should have retained their dark halos but observations suggest 
that globular clusters are not embedded in dark matter halos (Moore 1996; Heggie h Hut 1995). 

Another related population is the nine dwarf spheroidals in the Local Group with central 
velocity dispersions a <^ 10 km s _1 , including five below 7 km s _1 (e.g., Mateo 1998). In the 
hierarchical clustering scenario, the dark matter in a present halo was most probably divided at 
reionization among several progenitors which have since merged. The velocity dispersions of these 
progenitors were likely even lower than that of the final halo. Thus the dwarf galaxies could not 
have formed stars at high redshifts, and their formation presents an intriguing puzzle. There are two 
possible solutions to this puzzle, (i) the ionizing background dropped dramatically at low redshifts, 
allowing the dwarf galaxies to form at z ^ 1, or (ii) the measured stellar velocity dispersions of the 
dwarf galaxies are well below the velocity dispersions of their dark matter halos. 

Unlike globular clusters, the dwarf spheroidal galaxies are dark matter dominated. The dark 
halo of a present-day dwarf galaxy may have virialized at high redshifts but accreted its gas at 
low redshift from the IGM. However, for dark matter halos accumulating primordial gas, Kepner, 
Babul, & Spergel (1997) found that even if hi^z) declines as (1 + z) A below z = 3, only halos 
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with V c ^ 20 km s can form atomic hydrogen by z = 1, and V c ^ 25 km s is required to form 
molecular hydrogen. 

Alternatively, the dwarf dark halos could have accreted cold gas at low redshift from a larger 
host galaxy rather than from the IGM. As long as the dwarf halos join their host galaxy at a redshift 
much lower than their formation redshift, they will survive disruption due to their high densities. 
The subsequent accretion of gas could result from passages of the dwarf halos through the gaseous 
tidal tail of a merger event or through the disk of the parent galaxy. In this case, retainment of cold, 
dense, and possibly metal enriched gas against heating by the UV background requires a shallower 
potential well than accumulating warm gas from the IGM. Simulations of galaxy encounters (Barnes 
& Hernquist 1992; Elmegreen, Kaufman, & Thomasson 1993) have found that dwarf galaxies could 
form but with small amounts of dark matter. However, the initial conditions of these simulations 
assumed parent galaxies with a smooth dark matter distribution rather than clumpy halos with 
dense sub-halos inside them. Simulations by Klypin et al. (1999) suggest that galaxy halos may 
have large numbers of dark matter satellites, most of which have no associated stars. If true, this 
implies that the dwarf spheroidal galaxies might be explained even if only a small fraction of dwarf 
dark halos accreted gas and formed stars. 

A common origin for the Milky Way's dwarf satellites (and a number of halo globular clusters), 
as remnants of larger galaxies accreted by the Milky Way galaxy, has been suggested on independent 
grounds. These satellites appear to lie along two (e.g., Majewski 1994) or more (Lynden-Bell & 
Lynden-Bell 1995, Fusi-Pecci et al. 1995) polar great circles. The star formation history of the dwarf 
galaxies (e.g., Grebel 1998) constrains their merger history, and implies that the fragmentation 
responsible for their appearance must have occured early in order to be consistent with the variation 
in stellar populations among the supposed fragments (Unavane, Wyse, & Gilmore 1996; Olszewski 
1998). Observations of interacting galaxies (outside the Local Group) also suggest the formation 
of "tidal dwarf galaxies" (e.g., Due & Mirabel 1997). 

Finally, there exists the possibility that the measured velocity dispersion of stars in the dwarf 
spheroidals underestimates the velocity dispersion of their dark halos. Assuming that the stars are 
in equilibrium, their velocity dispersion could be lower than that of the halo if the mass profile 
is shallower than isothermal beyond the stellar core radius. As discussed in §2, halo profiles are 
thought to vary from being shallow in a central core to being steeper than isothermal at larger 
distances. The velocity dispersion and mass to light ratio of a dwarf spheroidal could also appear 
high if it is non-spherical or the stellar orbits are anisotropic. Some dwarf spheroidals may even 
not be dark matter dominated if they are tidally disrupted (e.g., Kroupa 1997). The observed 
properties of dwarf spheroidals require a central mass density of order O.IMq pc -3 (e.g., Mateo 
1998), which is ~ 7 x 10 5 times the present critical density. The stars therefore reside either in 
high-redshift halos or in the very central parts of low redshift halos. Detailed observations of the 
velocity dispersion profiles of these stars could be used to discriminate between these possibilities. 
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5. Conclusions 

We have shown that the photoionizing background radiation which filled the Universe during 
reionization likely boiled most of the virialized gas out of CDM halos at that time. The evaporation 
process probably lasted of order a Hubble time due to the gradual increase in the UV background as 
the HII regions around individual sources overlapped and percolated until the radiation field inside 
them grew up to its cosmic value - amounting to the full contribution of sources from the entire 
Hubble volume. The precise reionization history depends on the unknown star formation efficiency 
and the potential existence of mini-quasars in newly formed halos (Haiman & Loeb 1998a). 

The total fraction of the cosmic baryons which participate in the evaporation process depends 
on the reionization redshift, the ionizing intensity, and the cosmological parameters, but is not 
very sensitive to the precise gas and dark matter profiles of the halos. The central core of halos 
is typically shielded from the external ionizing radiation by the surrounding gas, but this core 
typically contains < 20% of the halo gas and has only a weak effect on the global behavior of the 
gas. We have found that halos are disrupted up to a circular velocity V c ~ 13 km s _1 for a shallow, 
quasar- like spectrum, or V c ~ 11 km s" 1 for a stellar spectrum, assuming the photoionizing sources 
build up a density of ionizing photons comparable to the mean cosmological density of baryons. 
At this photoionizing intensity, the value of the circular velocity threshold is nearly independent of 
redshift. The corresponding halo mass changes, however, from ~ 10 8 Mq at z = 5 to ~ 10 7 Mq at 
z = 20, assuming a shallow ionizing spectrum. 

Based on these findings, we expect that both globular clusters and Local Group dwarf galaxies 
with velocity dispersions <; 10 km s _1 formed at low redshift, most probably inside larger galaxies. 
The latter possibility has been suggested previously for the Milky Way's dwarf satellites based on 
their location along polar great circles. 
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David Spergel for useful discussions. We also thank Renyue Cen and Jordi Miralda-Escude for 
assistance with the reaction and cooling rates. RB acknowledges support from Institute Funds. 
This work was supported in part by the NASA NAG 5-7039 grant (for AL). 

APPENDIX A: Halo profile 

We follow the prescription of NFW for obtaining the density profiles of dark matter halos, but 
instead of adopting a constant overdensity of 200 we use the fitting formula of Bryan & Norman 
(1998) for the virial overdensity: 

A c = 18vr 2 + 82d - 39d 2 (10) 
for a flat Universe with a cosmological constant and 



A c = 18vr 2 + 60d - 32d 2 (11) 
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for an open Universe, where d = £l(z) — 1. Given Qq and we define 

am = n (i + ^) 3 (12) 

In equation @ c is determined for a given 8 C by the relation 

' ^ (13) 



3 ln(l + c) - c/(l + c) ' 
The characteristic density is given by 

5 C = C(/M,)(li^) 3 . (14) 

For a given halo of mass M, the collapse redshift z co u is defined as the time at which a mass M/2 
was first contained in progenitors more massive than some fraction / of M. This is computed using 
the extended Press-Schechter formalism (e.g. Lacey & Cole 1993). NFW find that / = 0.01 fits 
their z = simulation results best. Since we are interested in high redshifts when mergers are very 
frequent, we adopt the more natural / = 0.5 but also check the / = 0.01 case. [For example, the 
survival time of a z = 8, 5 x 10 7 Mq halo before it merges is ~ 30-40% of the age of the Universe 
at that redshift (Lacey & Cole 1993).] In both cases we adopt the normalization of NFW, which is 
C(0.5) = 2 x 10 4 and C(0.01) = 3 x 10 3 . 

APPENDIX B: Radiative Transfer 

We neglect atomic transitions of helium atoms in the radiative transfer calculation. We only 
consider halos for which ksT is well below the ionization energy of hydrogen, and so following 
Tajiri & Umemura (1998) we assume that recombinations to excited levels do not result in further 
ionizations. On the other hand, recombinations to the ground state result in the emission of ionizing 
photons all of which are in a narrow frequency band just above the Lyman limit frequency v = v^. 
We follow separately these emitted photons and the external incoming radiation. The external 
photons undergo absorption with an optical depth at the Lyman limit determined by 

dr 

— — = a m \y L )n m . (15) 
ds 

The emitted photons near vl are propagated by the equation of radiative transfer, 

= -cwi(y)nHlIv + r] u ■ (16) 

ds 

Assuming all emitted photons are just above v = ul, we can set aHii^) = OHiiyL) in this equation 
and propagate the total number flux of ionizing photons, 

F! = / ^dv . (17) 
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The emissivity term for this quantity is 

/ 



-T-dv = —a H in H iin e , (18) 
VL hv 4tt 



where olhi is the total recombination coefficient to all bound levels of hydrogen and uj is the fraction 
of recombinations to the ground state. In terms of Table 5.2 of Spitzer (1978), uj = (4>i — (p2)/4>i- 
We find that a convenient fitting formula up to 64,000 K, accurate to 2%, is (with T in K) 

uj = 0.205 - 0.0266 ln(T) + 0.0049 In 2 (T) . (19) 

When these photons are emitted they carry away the kinetic energy of the absorbed electron. 
When the photons are re-absorbed at some distance from where they were emitted, they heat the 
gas with this extra energy. Since ksT <C hvL we do not need to compute the exact frequency 
distribution of these photons. Instead we solve a single radiative transfer equation for the total flux 
of energy (above the ionization energy of hydrogen) in these photons, 

f°° I 

F 2 = / -^{hv-hv L )dv . (20) 

The emissivity term for radiative transfer of F2 is 

f 00 ^ 2.07 x IP' 11 Xl (p)- X2 (f3) _ 3 _ 1 

J ^\hv-hv L )dv= —q^ — n HII n e erg cm s sr , (21) 

where (3 = hui/kT, T is in K, and the functions %i and X2 are given in Table 6.2 of Spitzer (1978). 
We find a fitting formula up to 64, 000 K, accurate to 2% (with T in K): 

m - / 0-78 if T < 103 k 

Xi ( ) -X2{ ) - j _ Q 172 + 0.255 ln(T) - 0.0171 ln 2 (T) otherwise. ( ' 

From each point we integrate along all lines of sight to find t Vl , F\ and F 2 as a function of angle. 
Because of spherical symmetry, we do this only at each radius, and the angular dependence only 
involves 0, the angle relative to the radial direction. We then integrate to find the photoionization 
rate. For each atomic species, the rate is 

r4n poo T 

r 7i = / dn ^oi{v)dv s- 1 , (23) 

Jo Jvi hv 

where Vi and Oi(y) are the threshold frequency and cross section for photoionization of species 
i, given in Osterbrock [1989; see Eq. (2.31) for HI, Hel and Hell]. For the external photons 
the UV intensity is I v $e~ Tv , with the boundary intensity i^o = Iv L ,o{v l v L)^ a as before, and t v 
approximated as t 1/l (u/ul)^ 3 - Since ai(v) has the simple form of a sum of two power laws, the 
frequency integral in T 7 j can be done analytically, and only the angular integration is computed 
numerically (cf. the similar but simpler calculation of Tajiri & Umemura 1998). There is an 
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additional contribution to photoionization for HI only, from the emitted photons just above v = Vl, 
given by f ^d£l (Ji(i/£,)Fi. The photoheating rate per unit volume is n^i, where rij is the number 
density of species i and 



/•4tt roo J 

/ dD, / -^-(Ji{v)(hv — hvt)dv s _1 ergs s _1 . (24) 
Jo Jvi hv 



The rate for the external UV radiation is calculated for each atomic species similarly to the calcu- 
lation of I\yj. The emitted photons contribute to em an extra amount of Jq 77 dO G%(yi)F2. 



REFERENCES 

Abel, T., k Mo, H. J. 1998, ApJ, 494, L151 

Abel, T., Norman, M. L., k Madau, P. 1999, submitted to ApJL, |astro-ph/9812151 

Babul, A., k Rees, M. 1992, MNRAS, 253, 31 

Barnes, J. E., k Hernquist, L. 1992, Nature, 360, 715 

Bertoldi, F. 1989, ApJ, 346, 735 

Bertoldi, F., k McKee, C. F. 1990, ApJ, 354, 529 

Black, J. H. 1981, MNRAS, 197, 553 

Bond, J. R., Szalay, A. S., & Silk, J. 1988, ApJ, 324, 627 



Bryan, G. L., Machacek, M., Anninos, P., & Norman, M. L. 1998, to appear in ApJ, [astro- 



ph/980534C 



Bryan, G., & Norman, M. 1998, ApJ, 495, 80 

Burkert, A. 1995, ApJ, 447, L25 

Crotts, A. P. S., k Fang, Y. 1998, ApJ, 502, 16 

Due, P.-A., & Mirabel, I. F. 1997, Proceedings of IAU Symp. 186 (Kyoto), |astro-ph/97TT253 
Dekel, A., & Silk, J. 1986, ApJ, 303, 39 

Dey, A., Spinrad, H., Stern, D., Graham, J. R., & Chaffee, F. H. 1998, ApJ, 498, L93 
Donahue, M., k Shull, J. M. 1987, ApJ, 323, L13 

Elmegreen, B. G., Kaufman, M., k Thomasson, M. 1993, ApJ, 412, 90 
Efstathiou, G. 1992, MNRAS, 256, 43 



-20- 



Fall, S. M., & Rees, M. J. 1985, ApJ, 298, 18 

Ferland, G. J., Peterson, B. M., Home, K., Welsh, W. F., & Nahar, S. N. 1992, ApJ, 387, 95 
Franco, J., Ferrara, A., Roczyska, M., Tenorio-Tagle, G., & Cox, D. P. 1993, ApJ, 407, 100 
Fusi-Pecci, F., Ballazzini, M., Cacciari, C, & Ferraro, F. R. 1995, AJ, 100, 1664 
Garnavich, P. M., et al. 1998, ApJ, 509, 74 
Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 

Grebel, E. 1998, invited review to appear in IAU Symp. 192, "The Stellar Content of the Local 
Group" , |astro-ph/9812443| 

Haiman, Z., & Loeb, A. 1998a, ApJ, 503, 505 

— . 1998b, ApJ, in press, ^stro-ph/9807070| 

— . 1998c, invited contribution to Proc. of 9th Annual October Astrophysics 



Conference in Maryland, "After the Dark Ages: When Galaxies Were Young (the Universe 



at 2 < z < 5)", College Park, October 1998, |astro-pfi/9811395 



Haiman, Z., Rees, M., & Loeb, A. 1996, ApJ, 476, 458; erratum, ApJ, 484, 985 
Haiman, Z., Thoul, A. A., & Loeb, A. 1996, ApJ, 464, 523 
Haehnelt, M. G., & Steinmetz, M. 1998, MNRAS, 298, 21 

Heggie, D. C, & Hut, P. 1995, in IAU Symp. 174, Dynamical Evolution of Star Clusters- 
Confrontation of Theory and Observations, ed. P. Hut &: J. Makino (Dordrecht: Kluwer) 

Hu, E. M., Cowie, L. L., & McMahon, R. G. 1998, ApJ, 502, L99 

Katz, N., Weinberg, D. H., & Hernquist, L. 1996, ApJS, 105, 19 

Kepner, J. V., Babul, A., & Spergel, D. N. 1997, ApJ, 487, 61 

Klypin, A. A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, submitted to ApJ, |astro- 



ph/990124C 



Kravtsov, A. V., Klypin, A. A., Bullock, J. S., & Primack, J. R. 1998, ApJ, 502, 48 

Kroupa, P. 1997, NewA 2, 139 

Larson, R. B. 1974, MNRAS, 271, 676L 

Lacey, C. G., & Cole, S. M. 1993, MNRAS, 262, 627 

Lynden-Bell, D., & Lynden-Bell, R. M. 1995, MNRAS, 275, 429 



-21 - 



Majewski, S. R. 1994, ApJ, 431, L17 
Mateo, M. 1998, ARAA, 36, 435 

Miller, B. W., Whitmore, B. C, Schweizer, F., & Fall, S. M. 1997, AJ, 114, 2381 
Miralda-Escude, J. 1998, ApJ, 501, 15 



Miralda-Escude, J., Haehnelt, M., & Rees, M. J. 1998, submitted to ApJ, |astro-ph/98 12306 
Miralda-Escude, J., & Rees, M. J. 1998, ApJ, 497, 21 
Moore, B. 1996, ApJ, 461, L13 

Moore, B., Governato, F., Quinn, T., Stadel, J., & Lake, G. 1998, ApJ, 499, L5 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 (NFW) 
Navarro, J. F., & Steinmetz, M. 1997, ApJ, 478, 13 

Olszewski, E. W. 1998, in Galactic Halos: A UC Santa Cruz Workshop, ed. D. Zaritski (San 
Francisco: ASP) 

Peebles, P. J. E. 1993, Principles of Physical Cosmology (Princeton: Princeton Univ. Press) 
Peebles, P. J. E., & Dicke, R. H. 1968, ApJ, 154, 891 
Pen, U.-L. 1998, ApJ, 498, 60 

Petry, C. E., Impey, C. D., & Foltz, C. B. 1998, ApJ, 494, 60 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Quinn, T., Katz, N., & Efstathiou, G. 1996, MNRAS 278, L49 

Razoumov, A., & Scott, D. 1999, submitted to MNRAS, |astro-ph/9810425 

Rees, M. J. 1986, MNRAS, 218, 25 

Shapiro, P. R., & Kang, H. 1987, ApJ, 318, 32 

Shapiro, P. R., Raga, A. C, and Mellema, G. 1998, in Molecular Hydrogen in the Early Universe, 
eds. E. Corbelli, D. Galli, and F. Palla, Memorie Delia Societa Astronomica Italiana, 69, 
pp. 463- 469, |astro-ph/9804117 

Songaila, A., & Cowie, L. L. 1996, AJ, 112, 335 

Stecher, T. P., & Williams, D. A. 1967, ApJ, 149, L29 

Spitzer, L., Jr. 1978, Physical Processes in the Interstellar Medium (New York: Wiley) 



- 22 - 



Spitzer, L., Jr., & Hart, M. H. 1971, ApJ, 166, 483 
Tajiri, Y., & Umemura, M. 1998, ApJ, 502, 59 

Tytler, D. et al. 1995, in QSO Absorption Lines, ed. G. Meylan (Berlin: Springer) 

Thoul, A. A., & Weinberg, D. H. 1996, ApJ, 465, 608 

Unavane, M., Wyse, R. F. G.,& Gilmore, G. 1996, MNRAS, 278, 727 

Vader, J. P. 1986, ApJ, 305, 669 

Vader, J. R 1987, ApJ, 317, 128 

Verner, D. A., & Ferland, G. J. 1992, ApJS, 103, 467 

Voronov, G. S. 1997, ADNDT, 65, 1 

Weinberg, D. H., Hernquist, L., & Katz, N. 1997, ApJ, 477, 8 



This preprint was prepared with the AAS IAT^X macros v4.0. 



-23- 




Fig. 1. — Temperature T versus baryon overdensity A&. Dotted curve: tfj = t coo i, no radiation 
field. Solid curve: theat = tcooh with radiation. Dashed curves: tjj =net cooling/heating time, with 
radiation. The curves assume = 0.3, = 0.7, z = 8, a = 1.8 and rij/rib = 1. 
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Fig. 2. — Halo structure after reionization, for a halo mass M = 3 x 10 Mq at z = 8. Bottom 
plot: Baryon overdensity (A&) versus r/r v ; r . Middle plot: Neutral hydrogen fraction (uni/nii) 
versus r/r v ; r . Top plot: Ratio of thermal energy (TE) to potential energy (PE) versus r/r Y - ir . The 
dashed curves assume an optically thin halo, while the solid curves include radiative transfer and 
self-shielding. All the curves assume Qq = 0.3, J7a = 0.7, a = 1.8 and n^jn^ = 1. 
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Fig. 3. — Unbound gas fraction versus total halo mass. There are three pairs of curves, each 
consisting of a solid line (with radiative transfer) and a dashed line (without radiative transfer). 
From right to left, the three sets of curves correspond to a = 1.8, z = 8; a = 5, z = 8; and a = 1.8, 
z = 20. All the curves assume Q,q = 0.3, = 0.7, and n 7 /nfe = 1. 
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Fig. 4. — Total halo mass for which 50% of the gas is unbound versus reionization redshift. The 
solid line assumes a = 1.8, and the dotted line assumes a = 5, both for an S7o = 0.3, = 0.7 
cosmology. The other lines assume a = 1.8 but different cosmologies. The short-dashed line 
assumes CIq = 0.3, fl\ = and the long-dashed line assumes = 1- All assume n^/nf, = 1. 
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Fig. 5. — Circular halo velocity at which 50% of the gas is unbound, as a function of ionizing 
spectrum. Both panels show V c on the vertical axis. The left panel varies the spectral slope a for 
two values of the normalization, n-y/rif, = 1 (solid curve) and n 7 /nf, = 0.1 (dashed curve). The right 
panel varies the normalization for two spectral slopes, a = 1.8 (solid curve) and a = 5 (dashed 
curve). All curves assume an = 0.3, = 0.7 cosmology. 
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Fig. 6. — Fraction of the collapsed gas which evaporates from halos at reionization, versus the 
reionization redshift. The solid line assumes a = 1.8, and the dotted line assumes a = 5, both for 
an Qq = 0.3, = 0.7 cosmology. The other lines assume a = 1.8 but different cosmologies. The 
short-dashed line assumes f^o = 0.3, = and the long-dashed line assumes f2o = 1- All assume 
n^/rib = 1. 
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Fig. 7. — Total fraction of gas in the Universe which evaporates from halos at reionization, versus 
the reionization redshift. The solid line assumes a = 1.8, and the dotted line assumes a = 5, both 
for an Qq = 0.3, Oa = 0.7 cosmology. The other lines assume a = 1.8 but different cosmologies. 
The short-dashed line assumes fio = 0.3, Oa = and the long-dashed line assumes Qq = 1. All 
assume n^/nt, = 1. 



